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INTRODUCTION  | 

I ho  acceptance  of  optical  communications  for . in  operailona!  systems  has  been 
severely  huinpcied  by  our  inability  to  adequately  compcnsiite  lor  channel  effects  induced 
by  the  environment.  Consequently,  most  if  not  all  of  the  projected  system  gains  are  quick- 
ly nullified  when  rudimentary  measures  of  system  margin  are  added  to  the  link  budgets  to 
account  for  these  effects.  It  is.  therefore,  extremely  important  that  environmental  effects 
be  accurately  accounted  for.  and  systems  designed  to  best  exploit  these  channels  in  a most 
advantageous  manner.  TIic  most  difficult  channel  that  the  optical  communications  engine  ■ 
has  to  deal  with  is  the  multiple  scattering  channel.  Such  a channel  exists  in  propagating 
through  clouds,  fog,  water,  etc.  jref  1-3  ] . In  this  document  we  will  extend  the  model, 
which  has  been  independently  developed  by  Heggestad  |4]  and  Amush  |5]  for  multiple 
scattering  media,  and  apply  it  to  compute  the  effects  we  would  encounter  while  traversing 
u satellite  to  underwater  channel.  In  doing  so,  we  will  try  to  validate  the  use  and  interpre- 
tation of  the  model  by  applying  it  to  experimental  data. 

In  its  most  general  form,  the  problem  of  optical  communications  between  a satellite 
and  a submerged  platform  cun  be  described  as:  (I)  a problem  in  communications  from  u plat- 
form in  a nonscattering,  nondispersive  environment,  through  a random  surface,  and  into  a 
medium  with  a different  index,  which  is  multiple  scattering,  absorbing,  and  dispersive;  and, 
conversely,  (2)  a problem  in  communications  from  a platform  in  a multiple  scattering,  ab- 
sorbing, and  dispersive  medium,  through  a random  surface,  and  into  a medium  with  a dif- 
ferent index  which  is  nonscattering  and  nondispersive.  These  two  problems  are  nonrecipro- 
cal. Thus,  it  is  necessary  to  decompose  them  into  their  fundamental  elements  and  to  in- 
dividually identify  and  characterize  the  contributing  factors.  To  this  end,  this  discussion 
is  divided  into  four  parts.  The  first  part  (the  underwater  channel)  involves  the  actual  propa- 
gation effects  encountered  while  traversing  a multiple  scattering  medium.  The  three  system 
parameters  which  can  be  identified  are  attenuation,  beam  spreading,  and  apparent  source 
size.  These  in  turn  arc  related  to  the  absorption  coefficient,  the  scattering  coefficient,  and 
the  volume  scattering  function.  The  second  part  (the  air/sea  interface)  addresses  the  prob- 
lem of  transmitting  through  a random  surface  d aracterized  by  a slope  distribution.  The 
effect  on  scintillation  is  discussed  in  addition  t-  neam  pointing  and  beam  broadening.  The 
third  and  fourth  parts  address  link  calculations  from  the  satellite  platform  to  the  submerged 
platform  and  from  the  submerged  platform  to  the  satellite  platform,  respectively.  It  is 
estimated  that  tlie  model  presented  can  be  verified  to  within  several  dB  over  most  opera- 
tional scenarios  envisioned. 
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THE  UNDERWATER  CHANNEL 

Over  the  past  2 decades  there  has  been  an  interest  in  understanding  the  behavior  of 
light  while  propagating  through  water.  With  the  advent  of  the  laser  this  interest  intensified 
within  the  context  of  operational  equipment.  Although  there  have  been  numerous  measure- 
ments [6, 7 ] and  many  empirical  curves  have  been  derived  to  fit  the  data  (7  ] . the  latter 
are  of  limited  use  for  extrapolating  system  performance.  For  this  document,  we  will  use 
a model  which  has  been  developed  independently  by  two  separate  authors  (4, 5 1 . While 
this  model  is  derived  for  small-angle  forward  scattering,  it  appears  to  be  fairly  accurate  out 
to  ±4S”  provided  the  optical  thickness  is  neither  too  large  nor  too  small.  Fortunately,  the 
range  of  validity  is  within  the  operational  ranges  envisioned.  The  model  describes  the 
radiance  transfer  within  the  multiple-scattering  region.  This  region  is  characterized  by  three 
variables. 

1 . Absorption.  The  absorption  coefficient  of  the  medium,  a,  is  the  amount  of 
energy  absorbed  by  the  medium  per  unit  length  of  propagation.  This  loss  is  attenuation, 
and  goes  directly  into  heating  and  other  irreversible  processes. 

2.  Scattering  function  F(9).  Multiple  scattering  media  are  characterized  by  “mat- 
tering  centers  which,  in  the  case  of  ocean  water,  appear  to  be  both  from  plankton  and  from 
molecular  scattering.  The  volume  scattering  function  is  defined  as  the  secondary  radiation 
pattern  created  by  a plane  wave  traversing  a small  enougl.  volume  so  that  only  single  scat- 
tering occurs.  This  represents  the  average  scattering  distribution  of  all  the  scattering  centers. 
There  does  not  appear  to  be  a great  deal  of  variation  in  the  general  shape  of  F(0),  although 
the  average  width  does  change. 

3.  Scattering  coefficient.  If  we  normalize  F(0),  then 


1 

I 

1 

1 


JT 

s=2jt  J'  F(0)  sin  ddd  , (1) 

0 


and 


(2) 


is  the  normalized  version.  S is  the  scattering  coefficient  with  s~^  interpreted  as  the  average 
distance  between  scatterings.  Amush  assumed  a form  for  f(0)  as 


f(0)=-^e-^^  ;5s  10. 
2ir0 


(3) 


6.  Jerlov,  N.  G.,  "Optical  Oceanography,”  Elsevier,  1968 

7.  Duntley,  S.  Q.,  “Underwater  Lighting  by  Submerged  Lasers,”  Visibility  Laboratory,  SIO  ref  71-1 , 
June  197 1 
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lioggcstad.  on  tlie  other  hand,  defines  a modified  variance 


jr/2 

0-  = 7T  f 0-  sin0f(0)d0  ~ ^ , (4) 

0 

and  by  equating  0~  witli  !/6“.  the  two  models  are  identical.  This  is  true  for  large  5.  We 
will  use  the  notation  0~.  It  is  also  common  to  define  an  e.xtinction  coefficient  a.  defined 
by 


a = a + s. 


(5) 


Thus,  to  completely  characterize  the  environment,  it  would  be  necessary  to  have 
measuring  equipment  for  a.  s and  lYO).  An  alternate  procedure,  and  a less  desirable  one. 
would  be  to  measure  two  parameters  and  scale  the  measurements  to  the  third.  Although 
feasible,  this  would  assume  the  validity  of  a model  and  the  confidence  to  e.xtrapnlate  from 
, it.  With  these  parameters  in  mind,  we  will  now  present  the  model. 

For  convenience  we  assume  that  the  transmitted  beam  is  Gaussian  and  has  the  foim 


fo(0,_r)  = exp 

(itOqTq)- 


(6) 


that  is.  it  has  a Gaussian  distribution  in  both  its  spatial  cross  section  and  its  ray  direction. 
Next,  we  assume  the  geometry  in  figure  1 . where  the  source  is  at  (0.  0).  the  observer  is  at 
(_r,  z),  and  as  we  will  see.  the  apparent  source  is  at  (0.  zq).  In  terms  of  the  observation  point 
(r.  z)  we  have  as  the  transfer  in  intensity  f(0.  r) 


f(0,_r)  = exp 

(ttU^Rq)- 


% 


exp 

(7TU^R|)- 


(7) 


where 

1 + 2v  + +3  gy 

3(2  + V) 

=szO-(2  + V) 


r; 


= sz-'0- 


1 + 2V  + + 3CV 

2 + 3(«+V) 


T 


2 + 3(C  + V) 
_ 


(8) 
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Figure  1.  Model  geometry. 


(9) 


(10) 


Some  explanation  of  the  interpretation  of  equation  (7)  is  now  warranted.  First  notice  that 
if  we  had  a receiver  at  the  point  (r,  z)  and  added  equally  the  contributions  coming  from  all 
angles,  we  could  integrate  over  the  variables  8^,  6^  and  obtain  the  result* 


/ f<0-l)dn= 

(7tR“) 


(11) 


Thus,  the  total  energy  has  a distribution  which  is  Gaussian  in  the  z*plane,  is  centered  at 
r=0,  and  is  a result  of  scattering.  The  standard  deviation  of  this  spread  is  ((R ^ )/2)''^^  = 

|(sz^0“)/6  (2  + 3(8  + V)]  I If  in  addition  we  could  collect  all  the  scattered  radiation 
in  the  z-plane  (a  large  collector),  we  would  integrate  over  r and  obtain 


•Note;  If  the  ray  is  coming  from  the  directions^  = (0f,  6^,  then  the  receiver  is  pointed  in  the 
i-0)  = (-Oj.,  -0^)  direction.  Hence,  there  is  only  a sign  difference  between  the  two. 
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(12) 


f = y*  f(_r)dr  = e~^^'  , 

and  identify  this  as  an  irretrievable  loss  which  we  see  is  due  to  absorption. 

Now  suppose  we  observe  the  source  at  the  point  (r.  z)  as  a function  of  angle  (fig  2). 
Notice  is  the  unit  vector  in  the  direction  0^,  the  angle  out  of  the  r-z  plane,  and  Pj.  is 
the  unit  vector  representing  the  angular  tilt  up  from  the  z-a,xis  Oj.  in  a plane  described 
by£.  z.  Thus,  for  any  r,  the  maximum  always  occurs  in  the  r-z  plane  (0^  = 0)  at  a tilt  angle 
of  Alternatively,  for  a fixed  tilt  angle  0^.  the  maximum  occurs  when  the  receiver  is  off 
the  axis  a distance  rj^.  The  net  result  of  both  interpretations  is  that  the  source  appears  to 
be  located  at  the  point  zq  ( fig  1 ),  where 


Furthermore,  the  source  will  have  an  apparent  extent  (size)  in  diameter  (twice  the  standard 
deviation)  of 

I 3(  1 + V)  V 
or 

2(s0-z)'/-  (15) 

in  radians  (field  of  view).  Consequently,  any  system  should  account  for  the  spatial  filtering 
that  may  occur  when  optical  elements  are  used. 


•Si'" 


(14) 


Figure  2.  Source  as  a function  of  angle. 
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l inally.  we  can  identify 


rJ 


Mfi) 


as  the  ratio  of  the  initial  beam  cross  section  to  the  cross  section  at  (r,  i),  which  should  be 
much  less  than  one;  and  “ 


“ ■> 

szO"  U" 

r 


as  the  ratio  ot  the  initial  beam  spread  to  the  beam  spread  at  (r,  z),  which  should  also  be  much 
less  than  one.  Thus,  we  can  set  8 = V = 0 when  collimated  beams  are  used. 

Strictly  speaking,  the  model  used  here  is  only  valid  for  small-angle  forward  scatter. 
This  is  true  because,  in  the  derivation,  the  approximations  sin  and  cos  0~1  are  used. 
However,  these  approximations  are  only  off  by  10-20%  at  30-40°,  and  hence  the  model 
should  degrade  gracefully  at  larger  angles.  Some  modification  has  to  be  made,  however,  to 
use  this  at  large  angles.  This  is  due  to  the  fact  that  the  absorption  an^  scattering  paths  are 
longer  by  the  factor  zLse^O  - 1 1 at  the  angle  0.  This  can  easily  be  accounted  for  by  chang- 
ing z to  zsecO  = y//.~  +t'  wherever  z occurs.  Then  we  will  interpret  equation  ( 7)  to  be  the 
transfer  in  intensity  from  the  source  to  a sphere  of  radius  z.  With  the  latter  interpretation 
in  mind  we  will  now  show  the  justification  of  using  this  model  and  then  point  out  the  re- 
maining verification  needed. 

Consider  the  geometry  in  figure  3.  A collimated  source  emits  radiation  along  the 
/-axis.  Tile  medium  is  characterized  by  the  ratio  a/a.  Since  a = s + a,  s/a  = (a/a)  - 1.  The 
unit  of  length  is  N = a/,  extinction  lengths.  Thus,  N extinction  lengths  correspond  by  the 
relationship. 

N = az=^az=(-^)N3l, sorption  08) 

to  N/(a/a)  absorption  lengths  and.  since 


SPHERICAL  CAP  OF 
RADIUS;  AND  SOLID 
ANGLE  27rt1-COS);/] 


Figure  .3.  Geometry  used  by  Duntly  [7| . 
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Nj,cat  “ sz  = ((s+a)  - a]  z = az  - az  ■ 


to  N[(a/a)  - 1 !/(a/a)  scattering  lengths. 

With  this  geometry  and  parameterization,  Duntley  ' has  made  extensive  measure- 
ments of  the  power  collected  as  a function  of  N for  various  iues  of  Two  representative 
samples  are  shown  in  figure  4.  For  this  case,  equation  (1 1)  integrates  to 


y f(r)dr  = e"®^  Ll  - e * 


_ g-az  j _ g(2/3^s6^z 


1 -exp 


l(2/3)0 


As  pointed  out  by  Duntley  and  as  observed  in  equation  (12),  for  large  values  of  tp 
we  would  expect  the  relationship  e absorption  jq  jjqjj  However,  this  is  not  happening 
even  though  the  curve  saturates  for  ip  > 60*.  We  can  compensate  with  our  model  by  recall- 
ing that  the  effective  source  broadens  and  shifts  location  as  N increases.  Duntley  was  able 
to  observe  and  measure  the  former  phenomenon,  although  he  was  not  able  to  explain  it 
(fig  S).  In  the  actual  measurements  an  integrating  sphere  was  used  in  the  collecting  optics. 
This  has  a spatial  response  of  cos  0.  Consequently,  at  large  values  of  N one  would  expect 
to  start  to  observe  spatial  filtering  of  the  source.  This  is  precisely  what  we  see  for  N greater 
than  8.  This  was  corrected  for  as  follows.  It  was  assumed  that  the  f*iffereiice  between  the 

e"  absorption  and  the  p ^ 100°  curve  was  due  to  spatial  filtering.  Therefore,  at  every 
value  of  N this  difference  was  added  to  each  of  the  curves  (on  a log  scale)  and  replotted. 
The  model  was  then  calibrated  at  the  largest  value  of  N and  the  smallest  value  of  i/*,  where 
it  should  be  most  accurate,  and  O^was  calculated.  The  model  was  then  plotted  on  the  re- 
vised curves  (fig  6).  Although  the  agreement  is  not  perfect,  it  is  remarkably  close.  The 
validity  of  this  calibration  should  be  checked  at  some  point  by  comparing  the  results  of  an 
integrating  sphere  with  those  of  a hemispherical  coverage  lens  (fisheye). 
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Figure  6.  Revised  model  curves. 
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THE  AIR/SEA  INTERFACE 


In  this  section  a geometric  optics  model  is  developed  to  determine  the  effect  of 
surface  irregularities  on  beam  spreading,  pointing,  and  sci  itillation  in  traversing  the  boundary. 
To  do  so,  we  consider  the  following  model  of  an  element  of  the  surface  {fig  7).  A ray  of 
light,  y degrees  off  the  normal,  impinges  upon  the  surface  whose  local  slope  is  R degrees 
from  the  horizontal.  By  Snell’s  law 

n' sin  (7' + R)  = n sin  (7  + R)  (21) 

and 

7'  - sin-'  sin  (7  + R)j  - R , (22) 

which  is  valid  for  both  positive  and  negative  slopes.  7 and  7'  are  always  taken  with  respect 
to  the  true  vertical.  If  R is  a random  variable,  the  statistics  of  the  slope  must  also  be  fac- 
tored in.  We  do  this  in  the  following  manner.  Given  a sample  R from  the  set  of  possible 
slopes,  7'  is  well  defined.  That  is,  the  probability  of  7'  conditioned  upon  7 and  R is 

p(7'/7.  R)  = 6(7' -(sin"' sin  (7+R)^  -R)]  . (23) 

We  arrive  at  the  angular  distribution  0(7',  given  7,  by  averaging  over  the  variable  R.  Thus, 


00 

p(V/7)“  / dRpR(R)p(7'/7,R) 

•00 


00 

= f dR  pj^(R)6  I7' -(sin"' sin  (7+R)^  -R)l,  (24) 

—00 


where  pp(R)  is  the  probability  density  of  R. 

Rigorously,  this  is  merely  the  change  of  variables  in  the  density  Pr(R)  from  R to 
7'  by  use  of 


R = tan"' 


-7  sin  7 - sin  7' 
n 


> n 

cos  7 -~7  COS  7 
n 


= R(7'.7) 


(2.S) 
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NORMAL  TO 

LOCAL 

SURFACE 


VERTICAL 


lAY  PATH 


P(7'/7)=Pr  IR(7',7)I 

dy 


I 


With  7“  defined  ;is 


r>(7'/7)d7' 


= /|sin“‘  !^'~sin(7  + R)j  -r|  " pR(R)dR, 


the  rms  spread  becomes 


(28) 


—r:  —2  in 

A = (7'"-7'  ) '■ 


(29) 


There  are  some  practical  limitations  to  these  results  which  r<^quire  modification 
(fig  8).  A ray  of  light  with  zenith  angle  7 will  never  intercept  a wave  whose  slope  is  greater 
than  (ff/2)  -7  because  of  wave  obscuration.  However,  the  ray  will  still  penetrate  the  inter- 
face with  probability  one.  Consequently,  the  limits  of  integration  for  R are  set  at 
[ - (7r/2).  (r/2)  - 7I  and  the  density  Pr(R)  should  be  modified  to  that  of 


pr(R) 

ff/2-7 

/ PR(R)dR 
-*r/2 


p^R)  ; -ff/2<R<7r/2-7 
0 elsewhere . 


(30) 


The  results  in  equations  (26-29)  would  then  be  modified  by  replacing  pr(R)  with  p^(R). 

In  general,  the  results  presented  can  be  simplified  by  only  considering  those  values  off  + R) 
< 45“.  This  corresponds  to  the  major  operational  requirements  and  gives  good  engineering 
insight  into  the  behavior  of  a ray  going  through  the  air/sea  interface.  For  this  case 


Figures.  Modifications  of  model. 
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p(y'ly)  pr 


|l  -2,  I ^ varlRl  . 
n 


Notice  that  1 1 - n/n'l^  < 1 for  the  air/sea  interface  (index  of  water  = i.33,  index  of  air  = I), 
and  consequently  the  ray  spreading  is  appreciably  less  than_the  slope  spreading  of  the  ocean. 
In  addition,  the  surface  adds  the  contribution  of  (1  > n/n')R  to  the  normal  bending  due  to 
Snell's  law. 

In  all  cases,  the  model  used  represents  an  optical  beam  of  zero  cross  section  and  zero 
divergence.  The  exiting  beam  also  has  zero  cross  section  and  zero  divergence  but  is  being 
steered  by  the  roughness  of  the  surface.  If  this  surface  is  the  ocean,  piy’h)  represents  the 
average  time  history  of  the  beam  direction  with  p(y*ly)dy'  the  probability  that  it  is  point- 
ing within  a d^'  interval  of  the  y'  direction  at  any  instant  in  time.  This  apparent  beam  wan- 
der would  cause  severe  scintillation  in  an  operating  system.  As  the  cross  section  of  the  beam 
increases,  the  refracting  surface  can  no  longer  be  considered  locally  flat,  and  different  por- 
tions of  the  beam  are  refracted  at  different  angles.  Consequently,  it  would  be  possible  to 
average  out  the  beam  wander  in  the  direction  of  the  mean  7'  by  spreading  the  beam  over  a 
larger  portion  of  the  surface,  if  the  area  of  the  beam  is  A and  the  correlation  length  of  the 
surface  statistics  is  L,  then  there  are  approximately  A/(ffL/2)^  identically  distributed  inde- 
pendent paths  similar  to  diversity  paths.  If  we  further  assume  a depth  z such  that  the  beam 
cross  section  is  greater  than  A, 


(32) 


then  all  the  paths  will  overlap  at  the  receiver.  This  can  be  analysed  in  the  following  manner. 
First  we  notice  that  the  probability  of  having  the  beam  within  an  rms  deviation  about  the 
mean  is 

y + A 

f p(y'ly)dy'  (33) 

y-  A 


I 

I 


which  for  the  Gaussian  density  becomes  0.68.  Thus,  even  if  we  had  no  time  variations,  the 
beam  would  be  within  a deviation  of  the  mean  only  68%  of  the  time.  Now  suppose  we  pick 
N independent,  identically  distributed  paths  to  the  receiver  and  transmit  (1/N)^  of  the 
power  pj  in  each  path.  Since  the  paths  are  identically  distributed,  the  average  direction  of 
the  sum  is  still  y.  Now,  however,  the  variance  of  the  sum  becomes  A^/N,  or  a standard  de- 
viation for  the  sum  of  A/>/^  about  the  mean  7*.  If,  for  example,  we  set  N = 25,  and  assume 
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the  central  limit  is  approximately  valid,  the  probability  that  the  beam  is  within  is  now 
0.‘)99994.  Since  the  correlation  length  is  approximately  the  separation  between  independent 
spatial  Nyquist  samples,  we  see  that 


and  can  be  used  accordingly.  Furthermore,  it  can  be  shown  that  the  scintillation  will  reduce 
the  average  signal-to-noise  ratio  by  the  factor 


. (35) 

1 +a2/N 

A verification  of  these  results  and  the  relationship  in  equation  (32)  would  be  warranted. 

Finally,  we  can  interpret  the  function  piy'ly)  as  a beam-spreading  factor.  Thus,  if 
we  have  a propagating  beam  of  the  form  f(0,  r)  in  equation  (7),  or  f(7,  r),  then  the  output 
beam  after  traversing  the  surface  will  be 

f(7'.  r)=J  <iyp(y' ly)Hy,  r),  (36) 

or  an  average  over  all  input  ray  directions  weighted  by  the  relative  intensity.  Notice  that  we 
have  not  restricted  the  results  to  which  medium  corresponds  to  air  and  which  to  water. 

When  going  from  air  to  water,  set  n = I , n'  = 1 .33;  and  when  going  from  water  to  air,  set 
n « 1 .33  and  n'  = 1 . Then  the  computation  of  the  beam  moments  after  traversing  the  sur- 
face yields 


■y  = /7'f(7\r)d7' 
/f(7',  r)d7' 


var  l7'| 


-!(y'  -~y)~i\y\t)dy' 

/f(7'.  r)d7' 


(37) 


The  results  derived  in  this  section  were  performed  for  a one-dimensional  surface.  To 
extend  them  to  a two-dimensional  surface  is  straightforward  if  we  restrict  ourselves  to 
('artesian  coordinates.  The  variable  R would  then  become  the  pair  R.=  (x,  y),  and  the  one- 
dimensional results  would  carry  over  to  each  of  the  orthogonal  coordinates.  The  interpre- 
tation would  then  be  one  of  projecting  the  true  slope  distribution  onto  the  Cartesian  co- 
ordinate system.  Although  simple  in  theory,  the  actual  computations  are  difficult.  If  we 
use  the  linearization  implicit  in  equation  (31),  this  problem  is  greatly  simplified.  For  this 
reason  we  will  restrict  our  analysis  to  this  assumption.  To  make  the  calculations  for  large 
zenith  angles,  a more  rigorous  assessment  of  the  surface  geometry  must  be  performed  [81 . 


8.  Swennen,  J . P.  J . W.,  “Time-Average  Power-Density  Probability  Distribution  below  the  Ocean  Surface 
of  a Beam  of  Collimated  Optical  Radiation  Incident  on  the  Surface,"  JOSE  56,  p 224-229, 1966 
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SATELLITE  TO  SUBMERGED  PLATFORM 


The  computation  of  satellite  to  submerged  platform  power  budget  is  aided  by  a brief 
discussion  of  the  geometry.  It  is  assumed,  for  a variety  of  reasons,  that  we  will  project  a 
spot  on  the  ocean  approximately  1 mile  in  diameter.  Thus,  if  we  transmit  watts  of  radia' 
tion,  from  the  zenith  the  full  angle  of  the  beam  will  be  approximately  ( 1/22000)  radians 
^ 50  prad  = 10  s.  The  power  density  (intensity)  at  the  surface  will  be  approximately 

p 

L—  = 4.62  X lO""^  P(  W/m^.  (48) 

ir(830)“ 

If  the  surface  is  illuminated  at  an  angle  y from  the  zenith,  then  the  power  density  will  be 
Pt  c(  s 'y 

= 4,62  X 10-/  P,  cos  7 W/m^.  (39) 

ir(830)“ 

Now,  however,  the  circular  spot  has  elongated  into  an  ellipse  with  minor  axis  830  m and 
msyor  axis  (830/cos  yjm.  We  will  use  the  symmetry  along  the  mryor  axis  of  the  ellipse  to 
pick  a convenient  coordinate  system.  We  will  call  this  the  x-axis.  The  minor  axis  will  de- 
fine the  y-axis  of  the  coordinate  system  and  the  depth  of  the  ocean  will  constitute  the  z- 
axisd'ig  9). 

In  practice,  the  spot  will  have  a nonuniform  illumination.  We  will  account  for  this 
by  defining  a normalized  intensity  Kxq,  yg),  which  we  then  multiply  by  the  factor  In 
equation  (38).*  Notice  that  the  angle  7 is  always  measured  between  the  x-y  plane  and  a 
line  in  the  x-z  plane.  This  will  allow  us  to  use  the  second  form  in  equation  (7).  At  any  loca- 
tion (xq,  yg)  in  the  x-y  plane  an  elementary  surface  element  xq,  yQ  contributes  an  amount 
Kxq,  ygldxQdyQ.  A ray  with  this  intensity  passing  througli  the  air/sea  interface  yields  the 
value 

f(V;  ( X - xq),  (y  - yg))  = Kxq,  yg)  dxgdyg  p(7'/7)  (40) 

as  the  intensity  on  the  water  side  of  the  boundary.  At  this  point  we  will  consider  a func- 
tional form  for  p(V/7)  to  aid  in  the  computation.  From  experimental  results  [10]  it 


* To  be  correct,  the  transmission  coefficient  at  the  boundary  should  also  be  included  as  a function  of 
angle  [9). 

9.  Stratton,  J.  A.,  “Electromagnetic  Theory,"  McGraw  Hill,  1941 . 

10.  Cox,  C.,  and  Munk.W,,  Bulletin  Scrlpps  Inst  Oceanog  Univ  California j6  401, 19S6. 
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s 


, we 


find  at  a point  below  the  surface  that  the  contribution  from  the  intensity  I(xq,  Vq)  in  an 
area  dx^dyy  at  the  point  ( Xg,  yg»  0)  lo  Ihc  intensity  at  (x,  y,  z)  is 


AKx.  y,  z;y^7y)  = 


KxQ.yoldxodyo 


cxp  - 


a|z‘  + (xq  - x)“  + 'yg  -y)^I 


r + r 

2x  ^ 


R',2 


-201 


m 


where 


_ 3(1  + V*) 


2 + 3(fi'  + V')  n i 2 

Lv'z^+(xo-x)2  + (yo-y)2J 

R|^  = s|z“  +(Xq-x)“  + (yo“y)^)^'^^6^ 

U^^^slz^rtxQ- x)^+(yo-y)^J*/2  0^ 


(Xq-  X)(7x-7x)  (yp-  y)  ('Xy-'Ty)  ^ 

v/(Xq  - x)2  + (yg  - y)^  v/(Xg-X)2  + (yQ-y)^i 


^*x  *y 


’ (43) 


2 + 3(C'  + V')' 
3 


I + 2V'  + 68'  + 3E'V' 


2 + 3 (fi'  + V')  J 


(44) 


with 

2 

^ var  (R) 


S0“  (z“  +(Xg  - x)^  + (yg-y)“I  */“ 


^Xg^ 


=*  0 


s0‘  (z^  + (Xg  - X)“  +yg  -y)^)^/^ 


(45) 
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iP 

l 


( 

{, 


2 

X 


I 


y 


j-7x  + sin-' 


-7y  + sin“’ 


^Xq-x)  I - 
n/z-  + (Xq-X)^  i 

(yp  - y)  j “ 
N/z“  + (yQ-y)^j 


(z^  + (Xo-x)2) 


(z^+CyQ-y)-)  . 


(46) 


In  equation  (43)  we  have  introduced  ‘he  new  variables  and  7y.  The  former  is  the  angle 
measured  in  the  x-z  plane  while  the  latter  is  the  angle  perpendicular  to  the  x-z  plane.  This 
set  of  variables  results  from  a rotation  of  coordinates  of  the  variables  0^  and  0^  by  the 
transformation 


^yp-y) 


>/(yo-y)^  + (xq  - X)- 


Tx 


v^yp  - y)“  + (xp  - x)2 


0r  = 


(Xp  - x) 


(yp  - y) 


y)2  + (xp-x)2 


- y(yp  - y)^  + (xp  - x)- 


(47) 


The  necessity  for  this  rotation  arises  from  the  fact  that  Oj.  lies  in  the  plane  described  by  the 
points  (Xq,  yp),  (x.  y ) and  tjie  refracted  angle.  Consequently,  it  is  necessary  to  project  the 
angular  contributions  onto  a common  set  of  coordinates  before  integration.  Since  the  trans- 
formation is  unitary,  the  variables  7^  and  7y  are  still  normalized  to  one.  Finally,  wu  see  that 

00 

I(x,  y,  z;  7v  Ty)  = //  AI  (x,y,  z)dXQdyQ.  (48) 


It  is  evident  that,  even  with  the  simplifying  assumptions  used,  the  model  is  complicated. 
Therefore,  in  terms  of  an  experiment,  it  is  important  to  pick  a geometry  such  that  we  can 
make  further  simplifying  assumptions.  For  example,  an  experiment  using  the  sun  at  zenith 
would  have  7 = 0,  Kxp.  yp)  - constant,  x = y = 0.  Wc  could  r ’ 0 pick  a calm  day  so  that  we 
can  assume  = Rv  ~ “ 0-  (•  addition  we  collect  over  a sphere  with  a fi  Jicye  lens,  vve 

have  ^ 


1(0,  0.  z)  ^ 
Kxp,  yp,  7) 


00 

// 


1_ 

ffR'i  ~ 


exp 


I ^ - 

- jalz-  + x5  + 


ypi 


1/2 


(49) 


Setting  xg  = tq  cosp,  = rg  sinp  and  assuming 


. -1  ’‘0 
sin  ‘ 


v/z^  + X5  \/z^  + x^ 


and 


sin 


Vo 


^0 


n/z^  + Vq  \4“  + yo 


we  have 


alz^+r.)  '“+ TTTT/21 

(2/3)s02(z2^,2^3/2| 

g (2/3)sfl'^lz2  + r^)^^^ 


ro‘*’‘0 


(50) 


(Notice  that  tor  z = 0 the  approximation  1 ' = 0 does  not  hold.) 

The  power  coll«cted  at  depth  z will  be  merely  AKO,  0,  z),  where  A is  the  size  of  the 
collecting  aperture.  Consequently,  a measurement  of  1(0, 0,  z)/1(xq,  yg)  over  many  extinc- 
tion lengths  would  indicate  the  validity  of  extrapolating  the  model  to  great  depths. 

We  have  now  presented  three  separate  methods  for  computing  the  power  loss  to  a 
depth  2 when  the  source  is  at  the  zenith  and  no  other  effects  are  considered.  By  order  of 
expected  accuracy,  they  are  equation  (49),  equation  (50)  and  equation  (7),  when  the  beam 
radius  rg  is  considered  large.  In  figure  10  we  plot  equation  (SO)  as  a function  of  the  upper 
limit  of  integration.  Notice  that  in  all  cases  convergence  occurs  when  the  radius  is  approxi- 
mately z/3  for  02  i 0.0 1 . As  02  increases  from  0.01  to  0.1 1 , the  effective  surface  area 
increases  and  the  total  contribution  decreases.  A calculation  of  equation  (49)  was  also  made 
and  the  result  was  within  a few  percent  of  that  calculated  by  equation  (50)  for  0^  = 0.01 . 
Finally,  when  we  use  equation  (7)  with  rg  large,  it  can  easily  be  shown  that  1(0, 0,  z)/ 

Kxg,  yg)  is  merely  e"^^.  This  is  plotted  together  with  the  previous  results  in  figure  1 1 as 
a function  of  z.  At  300  m and  0-  = 0.01  the  difference  was  only  3 dB.  This  result  implies 
that  the  diffuse  reflection  coefficient  |6)  when  measured  at  the  zenith  is  approximately 
the  absorption  coefficient. 

In  a practi<’al  system,  one  will  encounter  background  noise  arising  from  the  sky  and 
the  sun.  When  this  occurs,  the  use  of  a 47t  steradian  collector  will  admit  an  unacceptable 
amount  of  noise  into  the  detector  circuitry.  For  these  cases  it  can  be  shown  that  to  opti- 
mize the  received  signal-to-noise  ratio,  a spatially  matched  filter  should  be  used.  Simply 
stated,  the  matched  filter  will  take  two  forms  depending  upon  whether  we  have  blue  sky 
oi  the  sun  (or  both).  To  eliminate  a source  such  as  the  sun,  the  filter  reduces  to  an  obscura- 
tion covering  the  field  of  view  subtended  by  the  sun  to  the  receiver.  For  an  extended  source, 
the  filter  takes  on  the  angular  distribution  subtended  by  the  source  to  the  receiver.  In  prac- 
tice this  reduces  to  an  obscuration  which  only  passes  that  portion  of  the  field  in  which  the 


24 


(0.0.Z) 


I (0. 0.  Z) 


1 X 10'^  I.  L,.  ,J  I I—  J L I 

10  20  30  40  50  60  7I 

DEPTH  IN  METERS 

Figure  1 1 . Equation  (7)  plotted  with  large. 
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m^or  portion  of  the  source  subtends.  Mathematically,  we  would  integrate  equation  (48) 
over  the  variables  7^  and  7y  with  the  integration  boundary  determined  by  the  receiver 
field  of  view.  Then  equation  (43)  could  be  rewritten  as 


AI(  X,  y,  z;  7'x.  7y) 


Kxq,  yp)  dxpdyp 
(trU^R',)2 


exp 


alz^  + (X  - xo)‘  + y - Vq)^  1 ^ 


*x  *y 


R'l^ 


1 


k' 


^0 


<XQ~x)o;n  \ 

-s/(xo-x)2  + (yQ-y)2y 

+ ,v.7  + _!5SL2l£in_t1 

''  y.xo:VHvo:7)^jJ 


(51) 


And,  if  we  assume  that  the  receiver  will  be  pointing  at  the  refracted  angle  (n/n')7.  we  can 
perform  the  integration  over  a finite  field  of  view  between  say  n/n'7  - A and  n/n'7  + A. 

If  we  perform  the  integration  over  a cone,  we  find  that  some  difficulty  would  arise  in  trying 
to  obtain  a closed-form  solution.  However,  by  referring  to  figure  12,  we  see  that  upper  and 
lower  bounds  can  easily  be  obtained  in  closed  form.  The  resultant  received  power  over  the 
finite  field  of  view  n can  be  obtained  by  using  combinations  of  the  function 


I(x,  y,z.  S2)=  f f 
XQ  Yq 


yo> 


irR 


>2 


exp 


(f^  + t^) 

alz2  + (x-xo)2  + (y-yo)2)l/2  + i|^^ 


G. 


(52) 
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Figure  12.  Integration  upper  and  lower  bonds 
In  closed  form. 


v/(x-xo)^  + (y-yorj 

1 r . /.  n\-=  1 . 


-^erf 
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Mx-XQr  + (y-yo) 
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1 1 r / n \ •=•  (yo~y^^m  1 

— erf  A -|1  -^  ) Ry  + m' 

“ ' " ' -v/(x-xo)^  + (y'yo>^ 
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For  small  fields  of  view,  Ci  can  be  replaced  by 


0 = JLexp--L  -(i-A) 

irV'J  Ul-  V V " / 


^XQ-X)0m 

/(xo-x)“  + {yo-y) 


<yo-y><^m 

4xq^  x)^  + (yo  - y>“ 


On  a final  note,  it  is  possible  to  obtain  an  estimate  of  the  pulse  spreading  by  referring  to 
figure  1 3 for  the  zenith  geometry.  If  the  primary  contributions  come  from  the  disc  with 
diameter  z,  then  the  maximum  path  difference  is 


p »"|(sec0  - I) : seed  * 1.1 15 


and  the  maximum  time  difference  is 

At=il.'p  = l^  (.1151  10->°z. 

c 2c 


(55) 


At  z “ 300  m,  AT  = 90  ns.  If  the  primary  contributions  come  from  twice  the  disk  diameter, 
AT  = 324  ns. 

Heggestad  has  computed  the  impulse  response  of  the  medium  from  which  he  evalu- 
ates the  delay  spread  as  the  I /e  point.  This  value  takes  the  form 


For  z = 300  m,  this  yields  193  ns 


(56) 


Figure  13.  Zenith  geometry  for  estimating 
pulse  spreading. 


SUBSURFACE  TO  SATELLITE  BUDGETS 


The  pari  of  the  system  most  difficult  to  model  has  been  the  siibsurfac;'  to  satellite  up- 
link. This  difficulty  can  best  be  understood  by  showing  why  the  two  are  not  reciprocal. 

On  the  downlink  a 1-miie  spot  projected  from  22000  miles  represents  an  antenna  gain  of 

4>L02p2£-<,6dB, 

On  the  uplink,  however,  if  w^o  through  one  scattering  length  of  water,  the  beam  solid 
angle  will  be  approximately  0-  or  an  antenna  gain  of  4ir/0^.  Since  10~“  - 10"^  the 
gain  is  only  21-31  dB.  The  gain  then  goes  down  as -10  log  in  scattering  lengths.  Be- 
cause of  the  paucity  In  gain  it  would  be  necessary  to  operate  the  system  closer  to  the  surface 
on  the  uplink  than  for  the  downlink,  if  the  scattered  radiation  as  described  by  the 
Heggestad-Arnush  approximation  is  used  exclusively. 

For  this  portion  of  the  link  it  is  necessary  to  investigate  the  radiation  in  greater 
detail.  To  do  so  it  is  helpful  to  use  the  normali/ed  version  of  the  Mutual  Coherence 
Function*  * *‘^*^Nspatial  covariance  function).  For  the  scattering  function  described  In 
equation  (.3)  this  becomes 

_ T 

-jr“ro“p“  j 

7(p)  = exp  - + sZ  — r ' * * 

■ v/i-KkoP)2fl2 

Notice  that  at  p = 0 this  is  normalized  to  unity,  and  we  assumed  a gaussian  source  with  an 
aperture  ecjual  to  7rr„“  focused  at  infinity.  In  normal  system  design  we  are  commonly 
interested  in  the  beamwidth  of  the  antenna  defined  at  the  3 dB  points.  Since  the  Mutual 
Coherence  Function  (MCF)  is  the  transform  of  the  angular  distribution  of  the  source  as 
seen  by  the  receiver,  we  can  equivalently  define  a ‘coherence  length’  p^  as  a comparable 
measure  of  antenna  eollimation.  Titus,  the  greater  the  coherence  length,  the  closer  the 
source  appears  to  approximate  an  impulse  in  angle  tpoint  source).  By  setting  the  MCF 
equal  to  e"  *’^^’^l-3dB)  and  solving  for  p^,.  we  can  investigate  the  behavior  or  the  radiation 
as  it  traverses  the  scattering  medium.  The  expression  for  p^.  becomes 


p.  = \ ) : Z< 

*■  \jrr,.  / s 


.38()(sZ)-.48j/0“] 
'2rr|^-;by3| 


which  is  shown  in  figure  14  as  a function  of  Z for  an  initial  divergence  of  10"^  radians,  and 
for  the  water  properties  defined  by  s = 0.6.  a = .04  and  0-  = .0 1 . Notice  that  for  a distance 
Z = .693/s  the  beam  propagates  as  it  would  in  vacuum,  and  the  correlation  length  increases 
as  the  beam  diverges.  However,  the  scattering  mechanisms  abruptly  take  hold  at  this  dis- 
tance and  the  coherence  length  decreases  dramatically  in  a very  short  distance,  and  rapidly 


* * Wells,  W.  H.,  "Loss  of  Resolution  in  Water  as  a Result  of  Multiple  Small-Angle  Scattering," 
JOSA,  vol  .S'),  no  (),  106‘).  p 1 10‘) 
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approaches  the  value 


sO-Z, 


(59) 


ticfinini;  the  Hceeestad-Arnush  approximation.  Since  this  behavior  is  dependent  upon  the 
scattering  properties  ol  the  water,  it  is  instructive  to  deline  the  albedo'"'  ce  as  the  ratio 
s/(s+a)  and  tile  extinction  length  as  N = aZ  = (a+s)Z  and  replot  figure  14  for  various  water 
parameters  in  figure  I 5.  Thus  we  see  that  we  rapidly  lose  the  gain  (or  imaging  capability) 
of  the  medium  as  we  traverse  a few  scattering  lengths,  which  can  vary  in  terms  of  the  extinc- 
tion length.  This,  however,  is  not  the  whole  story. 

If  we  observe  the  MCF.  eipiation  (57).  for  large  values  of  p.  we  observe  the  asymptotic 
value  of  e”^^.  From  Fourier  transform  theory  we  know  that  this  corresponds  to  a point 
source  which  relates  to  the  unscattered  portion  of  the  beam.  And.  while  the  power  associ- 
ated with  this  portion  of  the  beam  is  significantly  less  than  that  associated  with  the  scattered 
radiation  (lower  by  e“^^).  it  nevertheless  retains  the  full  gain  of  the  original  source.  Conse- 
quently. it  can  be  shown  that  for  the  uplink  geometry  a receiver  located  out  of  the  scattering 

media  at  a great  distance  from  the  source  will  always  collect  more  power  from  this  unseat- 

’ <1^1 

tered  component  than  from  the  scattered  component.  It  is  therefore  possible  to  consider 
the  uplink  radiation  as  composed  of  two  additive  (iaussian  terms  - the  first  retaining  all  the 
geometric  properties  of  the  radiated  source  but  attenuated  by  the  factor  and  the 

second  consisting  of  the  scattered  portion  of  the  radiation  as  considered  for  the  downlink. 
Father  of  these  terms  may  be  used  to  develop  a system,  but  the  resulting  systems  will  have 
vastly  differing  operating  scenarios  due  to  the  difference  in  coverage  and  pointing  require- 
ments. Consequently,  we  will  consider  the  general  problem  of  a Gaussian  beam  propagating 
up  through  the  air/sea  interface  and  determine  the  effects. 

We  will  first  make  the  computations  outlined  in  equation  (36)  for  the  Heggestad 
Arnush  approximation  and  then  show  how  the  Unscattered  result  follows. 

(This  component  is  actually  diminished  by  tlie  factor 


1 - exp 


+ ( x-nq'"  (y~yo^~ 


which  we  will  ignore.) 


To  compute  the  surface  irradiance  profile  upon  passing  through  the  interface  on 
the  uplink  it  is  only  necessary  to  insert  equation  (7)  into  equation  (36).  This  yields 


f(7-£)=  J*d7  p(7V7)  f(7.  r) 


(60) 


(we  assume  a collimated,  zero  cross  section  source  ((i!=V=0)).  F'inally.  to  detemiine  the 
angular  distribution  of  the  beam.  f(7'.  r)  is  integrated  over  the  surface  to  yield 

f(7')=  J*f(7'.  r)  tlr  = J*f(7,  r)  p (7/7)  d7dr.  (61) 


To  perform  the  integration  in  equation  ((i0).  several  factors  must  be  taken  into  account. 
First  recall  that  we  have  definetl  the  medium  containing  the  source  to  have  the  index  n and 
zenith  angle  7.  riius,  if  we  want  a more  figurative  description,  we  should  turn  the  coordi- 
nate system  in  figure  9 Lipside  down  to  yield  figure  16. 

Next  recall  that  we  projected  the  true  slope  statistics  of  the  surface  onto  the  x and  y 
coordinates.  However,  the  scattered  beam  has  circular  .symmetry  with  regard  to  the  angular 
divergence.  Consequently,  it  is  again  necessary  to  rotate  the  axis  of  the  angular  coordinates 
by  the  transformation  in  equation  (47),  which  will  allow  us  to  perform  the  integration  in 
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c(|iiation  (fiO)  in  C artesian  coordinates.  Taking  the  latter  remarks  into  account  and  again 
assuming  {Jaussian  slope  statistics,  we  obtain  the  function  f(7'.  r)  with 


fty.i)  =P{  j- 


dO^tlOy 


'Rl'I-r2ff|  1-2;  i 
L Ini 


varlRI 


R'r  V " Jj:?) 


/x2  + y2 


2 l-^l  var  [R1  L\  / 

0 "7)  *^y)  ]| 


where  we  have  assumed  that  the  source  is  located  at  (xq  ® 0,  Vq  = 0,  z)  and  the  surface  ir- 
radiance  profile  is  over  (x,  y,  z).  Consequently,  the  angular  distribution  of  the  emerging 
beam  is  obtained  by  integrating  ff7',£)  over  the  variables  (x,  y). 

As  the  surface  roughness  goerto  zero  in  equation  (62),  pCt’/t)  approaches  the  delta 
function  6(7'  - 7).  For  this  case  the  integration  can  be  performeffover  0^,  6^  to  yield 


f(7',  r)  = 


exp-  ixy/rr  x“  + y~  + ■ ^ ^ y ■ 


"-4-^ 

V " \/x-  + y-/  V " \/x-  + y-/  J 

We  can  also  perform  the  integration  in  equation  (62)  for  the  narrow-beam  case.  For 
this  case  we  can  extend  the  integration  from  -«>  to  <»,  yielding 

‘ I ~ y 1 I + x“  + y- 

/ L ^JjTrRi- 


c + f: 
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For  the  case  in  which  the  ocean  roughness  is  absent,  equation  (60),  we  sec  that  the 
exiting  beam  is  centered  around  the  linearized  Snell’s  angle,  tn^' >7^.  but  somewhat  steered 
toward  the  zenith.  Physicaily,  this  is  due  to  the  fact  that  the  scattered  ^^aths  whicii  lie 
closest  to  the  zenith  traverse  u shorter  distance  and  consequently  are  absorbed  less,  which 
.skews  tlie  beam. 

We  can  integrate  the  contribution  at  an  angle  2 over  the  entire  surface  as  indicated 
in  equation  (58).  This  results  in  the  function  ffy'),  which  is  the  angular  power  distribution. 
To  compute  an  uplink  budget,  otie  needs  to  integrate  ft^')  over  the  solid  angle  /subtended 
by  the  collecting  aperture,  At  a distance  R laigc  enough  for  the  function  ffV)  to  be  con- 
stant over  the  collecting  aperture,  this  solid  angle  is  merely  A/4)rR^.  Consequently,  the 
collected  power  is 

Pr  = f(y)-~  . (66? 

~ 4rrR2 

An  important  case  to  notice  is  when  the  beam  is  exiting  the  water  at  the  zenith  and  the  sur- 
face is  smooth.  Then  equation  (63)  integrated  over  the  hemisphere  yields  the  result  in 
equation  (4')),  multiplied  by  the  normalized  angular  distribution 


If  we  look  at  the  3-dli  contour,  the  beam  half  angle  is  (n/n’lU^  >/0  693.  and  the  effective 
gain  is 


The  link  loss  for  this  case  becomes  (7^  ~y'y~  where 
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Pt 


// 


iU, 


exp  - L|z- + X- + y~l  +.^5 — I 
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dxdy 


1 XQ  = yo  = 7x  = 0 
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(68) 


(69) 
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To  obtain  the  results  for  the  unscattered  beam,  we  would  use  equation  (6)  in  place 
of  equation  (7).  This,  however,  nierdy  requires  the  substitutions 


S = 0 


■7 


H't  -»  Pje"*^';  2'  = \/z^  x"  + 


(70) 


with  = 0.  With  these  substitutions,  equations  (62,  63.  64,  67,  68,  69)  become 

equations  (71  - 76),  respectively: 


»■(> 


. r)  = p / ^ - 


dd^dfly 


r?rZ'  0“)-2irll--2,rvar|Rl 

° n 


exn  ■ 


_L  2 j.  2 

+ X'^  + y'^ 


>2  ^ >2 

*y 

0^  (z“  + X 


- + y“)  ^5 


(71) 
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We  again  point  out  that  lor  large  zenith  angles,  the  iinearization  used  to  derive  these 
results  is  not  valid,  and  a more  detailed  analysis  is  required.  ^2 

Comparing  equations  (68)  and  (75)  and  letting  cos  X = 1 we  see  that  the  ratio 
of  the  gain  in  the  unscattered  to  scattered  beam  on  axis  is 


(X/rg)- 


where  we  have  assumed  that  the  beam  is  focused  at  infinity  (6q  = X/r^)  and  we  are  truns- 

“N. 

mitting  from  scattering  lengths.  Since  the  unscattered  beam  has  e **  times  the  power  of 
the  scattered  beam,  we  see  that  the  inequality 


--Ns 

N,6-e 


(X/ro)- 


determines  which  portion  dominates. 
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DISCUSSION 


In  this  document  we  have  developed  models  for  use  in  evaluating  the  performance  of 
the  duplex  subsurface  to  above-surface  optical  communications  systems.  We  will  now  brief- 
ly discuss  the  limitations  of  the  models  and  the  areas  of  applicability.  We  will  also  point  out 
relevant  areas  of  future  work.  The  first  aspect  of  the  model  is  the  estimate  of  underwater 
propagation.  The  Heggestad-Arnush  model  used  appears  to  have  all  the  attributes  necessary 
for  accurate  predictions.  Although  the  model  has  been  calibrated  to  existing  data,  an  inde- 
pendent verification  is  warranted.  The  results  of  such  an  effort  would  determine  whether 
or  not  further  refinement  is  necessary.  The  major  implication  of  this  propagation  model  is 
a clear  distinction  between  the  contributions  of  absorption  and  scattering  to  the  extinction 
coefficient.  Couched  in  system  terminology,  the  model  states  that  if  the  size  of  the  beam 
on  the  surface  is  comparable  to  or  greater  than  the  depth  from  which  it  is  to  be  viewed, 
and  if  the  field  of  view  at  the  receiver  can  be  made  large,  then  the  only  loss  Is  from  absorp- 
tion. Since  the  extinction  coefficient  is  usually  two  or  more  times  greater  than  the  absorp- 
tion coefficient,  the  depth  prediction  under  these  conditions  would  also  be  two  or  more 
times  greater  than  predicted  from  the  extinction  coefficient  alone.  The  field  of  view  en- 
countered is  also  seen  to  be  proportional  to  the  square  root  of  the  scattering  coefficient 
and  the  depth. 

The  second  aspect  of  the  model  relates  to  the  effects  of  the  surface.  In  order  to 
obtain  usable  results,  a linearization  of  Snell's  law  was  employed  which  should  be  reason- 
ably accurate  for  zenith  angles  out  to  45°.  The  important  result  is  that  at  sufficient  depth 
the  effects  of  a random  surface  would  be  negligible.  The  major  concern  should  be  blockage 
of  light  due  to  foam,  etc.  The  basis  of  this  conclusion  stems  from  the  prediction  that  the 
rnis  beam  spreading  will  be  proportional  to  II  - (n/n'l  times  the  rms  slope  distribution  of 
the  surface.  This  would  imply  that  a maximum  of  5°  or  so  is  all  that  would  ever  be  expect- 
ed. The  major  impact  would  seem  to  be  on  the  uplink,  where  beam  steering  would  occur. 

The  most  difficult  part  of  this  communication  system  appears  to  be  the  uplink. 
Because  of  the  nonreciprocal  nature  of  the  duplex  system,  the  unscattered  portion  of  the 
beam  provides  the  greater  potential  for  power  transfer.  However,  the  power  in  this  portion 
of  the  beam  is  greatly  diminished  over  the  scattered  term,  yet  retains  its  high  directionality. 
The  diminished  power  implies  a depth  reduction  of  a/(a+s),  In  addition  the  spot  size  on  the 
ocean  suilace  will  not  encompass  enough  area  to  average  out  the  dynamic  effects  of  the  wave 
motion.  Consequently,  measures  will  have  to  be  taken  to  compensate  for  this  wave  motion 
in  an  active  and  dynamic  manner  This  implies  a form  of  image  inhancement  of  the  down- 
link beam  so  as  to  track  tlie  unscattered  component.  This  is  an  area  where  future  work  can 
he  directed  and  efforts  are  already  underway. 
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